The Fast Wandering of Slow Birds 
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I study a single "slow" bird moving with a flock of birds of a different, and faster (or slower) 
species. I find that if the speed Vs of the "slow" bird equals a certain "magic" speed 112, different 
from the speed of the flock, it will randomly wander transverse to the mean direction of flock motion 
far faster than the other birds will: its RMS transverse displacement will grow with time t like , 
in contrast to for the other birds. If «s 7^ «2, the RMS displacement of the "slow" bird crosses 
over from te to fs scaling at a time tc that scales according to tc oc \vs — Vc\~'^ ■ 



Flocking [H - the collective, coherent motion of large 
numbers of organisms - is one of the most familiar and 
ubiquitous biological phenomena. This phenomenon also 
spans an enormous range of length scales: from kilome- 
ters (herds of wildebeest) to microns (e^., the microor- 
ganism Dictyostelium discoideum 0, y, 

Following the recognition that this phenomenon is 
a dynamical version of ferromagnetic ordering, a phe- 
nomenological continuum theory 0, 0, S, HI was de- 
veloped, which showed, among other things, that two- 
dimensional flocks did not obey the Mermin- Wagner the- 
orem 10]: that is, they are capable of spontaneously 
breaking a continuous symmetry (rotation invariance) by 
developing long-ranged order, even in flocks with only 
short ranged interactions. Such order is impossible (in 
spatial dimensions d = 2) for equilibrium systems [10|. 

This fundamental difference between flocks and equi- 
librium systems arises from novel "anomalous" scal- 
ings of fluctuations in flocks. Many of the quantitative 
predictions of the phenomenological continuum theory 
about these scalings have been confirmed in numerical 
simulations [7I [ill . In particular, the lateral wander- 
ing of a bird in a flock of identical birds was found to 
grow superdiffusively in d = 2, scaling like i^0jll2|- 

In this paper, I consider a single "slow" bird, different 
from the other birds in the flock, which attempts to follow 
the other birds, but flies at a speed 7^ wo, where vq 
is the average speed of the flock. I find that if Vs = V2, 
where V2 vqis a "magic" speed characteristic of the 
flock, the "slow" bird will wander laterally much faster 
even than the superdiffusive wandering of the others. 

Most of this paper deals with the two-dimensional case, 
for which the anomalous scaling exponents are known. 
In three dimensions, the scaling exponents are unknown, 
but a plausible, albeit unproven, conjecture for them ex- 
ists. At the end of the paper, I will simply quote the 
results that apply for the three dimensional case if this 
conjecture is true. I will, in addition, quote results for 
four and more dimensions, where a linearized version of 
the continuum theory is adequate [l4|. 

My result applies to any "ordered "flock: i.e., one in 
which the flocking organisms spontaneously pick a direc- 



tion to move together via purely short-ranged interac- 
tions that make neighbors tend to follow each other, but 
which do not pick out any a priori preferred direction for 
this motion. That is, the flocking spontaneously breaks 
rotation invariance, as equilibrium ferromagnetism does. 

One specific realization of such a system is one obeying 
the Vicsek algorithm , in which birds attempt to follow 
their near neighbors, but make mistakes. This algorithm 
has a parameter range in which the system "flocks" (i.e., 
has a non-zero average velocity: < v >^ 0, where the 
average is over all the birds in a very large flock (i.e., in 
the thermodynamic limit)). In this "ordered" parameter 
range, the continuum theory 0, 0, S, [§| applies. 

Now let us add to this ordered flock a single "slow" 
bird[l^ which tries to follow the other birds in the flock 
via some short-ranged interaction, but which moves at 
a different speed Vs than the other birds in the ffock. 
For example, one could imagine (or simulate) a Vicsek 
algorithm [5| in which the "slow" bird obeys the same 
algorithm as the others, except that its speed of motion 
between time steps differs from those of the other birds. 

Consider now the mean squared lateral displacement 
of the slow bird: 



(1) 



perpendicular to the mean direction of motion of the 
flock. (Here and throughout this paper, _L and || denote 
directions perpendicular to, and along, the direction of 
mean flock motion, respectively.) 

Since the mean x±— position of the slow bird obeys 



xj:it)=xjriO)+ / v^{t)dt 



(2) 



where v-^{t) is the _L velocity of the slow bird at time t, 
w'^{t; Vs) is given by 



dt' / dt"{v^{t')v^{t")) 



(3) 



Now I need to relate the perpendicular velocity of the 
slow bird to the position and time dependent velocity 
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field v± (r, t) of the flock. Since the slow bird moves in 
the same direction[l^ as the other birds in the flock, the 
required perpendicular velocity is given by: 



(4) 



where vq is the mean speed of the flock (in the sense that, 
< v{r,t) >= voX\\, with i|| the unit vector along the 
direction of mean flock motion), and rs(t) is the position 



of the slow bird at time t. This is given by 



•fsit) — fsiO) + Vstxw + Sxl{t)x^\ + Sx^{t)x± 



(5) 



where Sxl{t) and Sx^{t) are the deviations of the slow 
bird from uniform motion at speed Vs along x^^. 

Using (|4]) and ([5]), I see that the desired single bird 
autocorrelation function in ([3]) is 



vo 



{v^ {rs{0) + {vsf + fell {t')) i[| + 5xi_ {t')x^,t') (f,(0) + {v,t" + fe|| {t")) i|| + Sx^_ {t")x^,t")) 



(^^^ Ci_ {x^ it') - x^ it") , Vs\t' - t"\ + fell it') - fell it") , t' - t") 



(6) 



where Cj_ (f, t) is the real space and time perpendicular motion Vstxu that their effect on the velocity-velocity au- 



velocity field auto-correlation function. 



tocorrelation in ((25|) is negligible. This leaves the task of 



I'll assume (and have verified a posteriori) that both evaluating C± [R± — 0,R\\ — Vst,t). 
(5x|| and dx±_ are small enough compared to the average Expressing C'j_ in terms of its Fourier transform gives 



aSt) ^ {v^ it') it")) - (^) ' {R^ = 0, i?|| = vjt) = (^) 7 ^^^^ II (,-, .) , (7) 

I 

where I've defined 6t = t' — t", and the spatio-temporally relation function C± (g, uj) is predicted by the continuum 
perpendicular Fourier-transformed velocity-velocity cor- theory Q to be: 



C±iq,u;) 



A (g; - V2q\\y 

iuj - c+ ie^) q)^ (a; - c_ (%) q)^ + {lo (Fl (g) + Fp iq)) - q\\ (i;2Fl (q) + iV p iq)))' 

I 



(8) 



where A is a phenomenological parameter related to the 
RMS magnitude of the directional errors the birds make 
in trying to follow each other, and the two direction- 
dependent sound speeds c± iO^) are given by 



C± {^^) = I I cos iOq) ± C2 (%) , (9) 



with 



V± ie^) = ± ( ^— ^ ) cos (%) + C2 (%) , (10) 



and 



C2 (0f) ^ 



_ Jil-V2) C0S2 (%) 



c2sin2(%) , (11) 



where 9^ is 



the angle between q and the direction of 
flock motion (i. e., the a;|| axis), and 7 and V2 are addi- 
tional phenomenological parameters with the dimensions 
of speed. Neither of them is generically equal to the mean 
speed Vo of the flock. The speed V2 will prove (as I'll show 
below) to be the "magic" speed that leads, if Vs = V2, to 
faster wandering of the "slow" bird; hence, the fact that 
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V2 7^ is crucial to making the behavior of the "slow" 
bird distinct from that of the other birds in the flock. 

In addition, the wavevector dependent longitudinal, 
and p dampings p are given by: 



Tl {q) = DL{q)ql 



(12) 




^ ' A ^ y A ) 
' A ^ y A ) 



(14) 



(13) 



where the renormalized, wavevector ((7)-dependent diffu- 
sion coefficient DL{q) obeys the scaling law: 



here A is an ultraviolet cutoff of order an inverse inter- 
bird distance. 



Performing the integral over u by complex contour 
techniques, and using eqn. ^ for C±{q,Lu), gives 



C{t) 



^1 A 



(Pq 



q^exp 



Vsqii 



r{0^)q)t- 



-{q)t^ {qCaidq) - "2911 



q'c2{e^)v_„{6^)T„{q) 



r 



(15) 



where the summation index a = 
modes with speeds c± [9^) ( Eqn 



T.{q) = Tl 



2C2 



± labels the two sound 
, (EI), and I've defined 



k) 



D^{q)ql) 



2C2 (% 
Vg (%) 
2C2 {9g-) 



D 



'P9|| 



2C2 



(16) 



Because of the anisotropy of the damping Fg., this in- 



tegral is dominated, as t — > oo, by g|| 
hence, 9^ — > 0, and I can replace Va 



{q±/A)^ A:^ qr, 
1^) and Cc{9^) ev- 



erywhere they appear with Va{9^ = 0) and c^{9^ = 0). I 
can likewise replace I?] = 9 with \q\\\. 

Doing this, we easily see that the scaling form eqn. ((T4|) 
for DL^q) implies a similar scaling form for Tc,{q)'. 



Ta{q) 




(17) 



Putting 9q- — in eqn. Q, I see immediately that 
one of the two c„{9^ = 0)'s is = W2, while the other 
= 7. Clearly, only the latter value of a survives, due to 
the Caq — V2q\\ factor in equation (jl5p . Keeping only the 
appropriate term then gives: 



C{t) = 2 — A 

"0 



dq^dqw e(^(^=-"2)9i|-r.(<r))t 



(2^)2 



,(18) 



where the value of a is that for which Co-( 



Now clearly, if Vg ^ "2, the first term in the expo- 
nential is much larger, at small g, than the ^^{q) term. 
Since the former scales as while the latter scales as 

q^, in the dominant regime q\\ ^ q^[. Hence, I can drop 
the latter term in the exponential, since it is small q's 
which dominate the integral at large times (as we'll see 
in a moment). Dropping this term, and using the scal- 
ing law p?|) for T^{q), it is easy to see that rescaling the 
components gy and q± of q according to 

q\\ = —, (l± = --r (19) 

t i 3 

pulls all of the time dependence out of the integral (Pq, 
giving the scaling law 

C(0 cxi-i (20) 

for C{t). Inserting this into ^ and ^ into ((3]) gives 

w'^{t;vs) (xt^ (21) 

a known 0] result for birds in flocks of identical birds. 

However, for slow birds, we now see that the results 
change dramatically if Vs ~ V2, since the {vg — V2)q\\t 
term in the exponential in eqn. (fT8|) vanishes. As a re- 
sult, I must keep the F^ term in the exponential in that 
equation. Doing so, and again using eqn. (fT7|) . I see that 
the rescaling to a new vector variable Q via 

1\\ ^ ^ (22) 



q± 



(23) 
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pulls all of the time dependence out of the integral when 
Vs = V2, and gives the scaling law 



C'it) = 



(24) 



where A is an unimportant constant. Note that this de- 
cay is much slower than that fean. (j20p )for a bird moving 
at a speed other than the magic speed Vg — V2- Note 
further that since, in general, the mean speed vq of the 
flock differs from V2 (i.e., vq ^ V2) the identical birds 
that make up the bulk of the flock will exhibit the more 
rapidly decaying correlations of eqn. (|20p . It is only a bird 
that differs from the generic members of the flock, specif- 
ically by flying at a speed Vs ^ fo, that has a chance to 
achieve Vs — V2, and the much slower decay eqn. (|24p . 
This slower decay of velocity correlations leads immedi- 
ately to faster lateral wandering. This can be seen by 
inserting into ([3]), which gives 



w'{t;vs = V2) = A / dt' / dt"\t'~t"\ 



9A 



(25) 



Since ^ ts ast^oo, this result eqn. p5|) confirms 
the fundamental conclusion of this paper: birds moving 
at the magic speed V2 wander laterally much faster than 
birds moving at any other speed. In particular, they 
wander much faster than the "normal" birds in the flock. 

What happens to a bird moving at a speed Vs close 
to, but not equal to, W2? For sufficiently short times, it 
should be impossible to tell that this bird was not moving 
at exactly the magic speed V2, and would therefore expect 
to recover the w'^ cx i 3 scaling I just found for that case. 
For longer times, though, one would expect the difference 
between the speed of the slow bird and the magic 
speed V2 to become apparent, leading to a crossover to 
the cx i 3 scaling found in that case. 

This proves to be precisely the case. To see this, note 
that the rescaling eqns. ((22)) . (|23p imply that the integral 
eqn (fT8|) for C{t; Vg) is dominated by parallel wavevectors 



^dom 



(X t 



(26) 



Comparing the Svq^^ term in the exponent of eqn (jlSp 
with the Tl term, where Sv = Vs — V2, and using the 
scaling law p?]) for and (|26p for the dominant q, I 
see that (5uqf?°™ ~ ^L{qdom) at a time tc which scales 



according to Svtc 



^ (x {tc ^)^, which implies 
tc{6v) oc Sv^'^ 



(27) 



For smaller times {t ^ tc), the F^ term in eqn. (fTSl) 
dominates the Svq\\ term, and the w'^{t) oc fs law found 
in eqn ((25l) holds. In the opposite limit {t ^ tc), the 
scaling w'^{t) cx ^3 holds. 

To summarize, for flocks in c? = 2 spatial dimensions. 



v?{t;vs 



cx 



t3^ t<t: tc{Sv) 
ti, t:s> tc{Sv) 



(28) 



with tciSv) given by ((27)) . 

The above results all hold in spatial dimension d — 2. 
In d ~ 3, unfortunately, the exponents characterizing 
scaling laws like eans. p^ . and (fT7|) are unknown. Fur- 
thermore, other quantities like A, and Dp could 
potentially also become strongly wavevector dependent 
in d = 3, as is in c? = 2. However, if it proves 
that only a particular one of the four potentially rele- 
vant non-linearities in the phenomenological equations 
of motion0,0,[l,[§| is actually relevant in d = 3 (namely, 
the "convective" one), which is certainly consistent with 
all calculations that have been done up to now, though 
by no means definitively established by them, then these 
complications disappear, and scaling laws analogous to 
eqns.(fT4|). and (fT7|) apply, but with different, but known, 
exponents (these are called the "canonical" exponents in 
[1,01). Assuming this to be the case, and repeating the 
analysis I've just done for the d = 2 case, I find, in d = 3, 



w'^{t;vs 



cx 



ti, t<$: tc{5v) 
t, t > tc{Sv) 



(29) 



with tc{Sv) again given by (|27|) . Note that the t 3> tc{Sv) 
behavior {w'^{t;Vs) cx t) is a simple "drunkard's walk". 

In the (potentially simulate-able) case of spatial dimen- 
sions d > 4, I find simple "drunkard's walk" wandering 
{w'^{t;Vs) cx t) independent of Vg- 
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